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Abstract. We study the stochastic Ising model on finite graphs with n vertices and 
bounded degree and analyze the effect of boundary conditions on the mixing time. 
We show that for all low enough temperatures, the spectral gap of the dynamics 
with (+)-boundary condition on a class of nonamenable graphs, is strictly positive 
uniformly in n. This implies that the mixing time grows at most linearly in n. 
The class of graphs we consider includes hyperbolic graphs with sufficiently high 
degree, where the best upper bound on the mixing time of the free boundary 
dynamics is polynomial in n, with exponent growing with the inverse temperature. 
In addition, we construct a graph in this class, for which the mixing time in the 
free boundary case is exponentially large in n. This provides a first example where 
the mixing time jumps from exponential to linear in n while passing from free to 
(+)-boundary condition. These results extend the analysis of Martinelli, Sinclair 
and Weitz to a wider class of nonamenable graphs. 



1. INTRODUCTION 

The goal of this paper is to analyze the effect of boundary conditions on the 
Glauber dynamics for the Ising model on nonamenable graphs. We will focus on 
a particular class of graphs which includes, among others, hyperbolic graphs with 
sufficiently high degree. Before discussing the motivation and the formulation of 
the results we shall give some necessary definitions. 

Given a finite graph G = (V,E), we consider spin configurations a = {cr x }xev 
which consist of an assignment of ±1 -values to each vertex of V. In the Ising model 
the probability of finding the system in a configuration o e {±l} y = &g is given 
by the Gibbs measure 




x€V 



where Z G is a normalizing constant, and (3 and h are parameters of the model 
corresponding, respectively, to the inverse temperature and to the external field. 
Boundary conditions can also be taken into account by fixing the spin values at 
some specified boundary vertices of G. The term free boundary is used to indicate 
that no boundary is specified. 

The Glauber dynamics for the Ising model on G is a (discrete or continuous time) 
Markov chain on the set of spin configurations Q G , reversible with respect to the 
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Gibbs measure \iq. The corresponding generator is given by 



where a x is the configuration obtained from a by a spin flip at the vertex x, and 
c x (a) is the jump rate from o to a x . 

Beyond of being the basis of Markov chain Monte Carlo algorithms, the Glauber 
dynamics provides a plausible model for the evolution of the underlying physical 
system toward the equilibrium. In both contexts, a central question is to determine 
the mixing time, i.e. the number of steps until the dynamics is close to its stationary 
measure. 

In the past decades a lot of efforts have been devoted to the study of the dy- 
namics for the classical Ising model, namely when G = G n is a cube of size n in 
the finite-dimensional lattice Z d , and a remarkable connection between the equi- 
librium and the dynamical phenomena has been pointed out. As an example, on 
finite n-vertex cubes with free boundary in Z d , when h = and j3 is smaller than 
the critical value /3 C (one-phase region), the mixing time is of order logn, while for 
(3 > (3 C (phase coexistence region) it is exp(n (rf-1 )/ d ) f 11511 12511251125m . 

More recently, an increasing attention has been devoted to the study of spin 
systems on graphs other than regular lattices. Among the various motivations 
which are beyond this new surge of interest, we stress that many new phenomena 
only appear when one considers graphs different from the Euclidean lattices, thus 
revealing the presence of an interplay between the geometry of the graph and the 
behavior of statistical systems. 

Here we are interested in the problem of the influence of boundary conditions 
on the mixing time. It has been conjectured that in the presence of (+)-boundary 
condition on regular boxes of the lattice Z d , the mixing time should remain at most 
polynomial in n for all temperatures rather than exp(n ( - d ~ 1 ^ d ) |j9t] . But even if some 
results supporting this conjecture have been achieved [53, a formal proof for the 
dynamics on the lattice is still missing. 

However a different scenario can appear if one replaces the classical lattice struc- 
ture with different graphs. The first rigorous result along this direction, has been 
obtained recently by Martinelli, Sinclair and Weitz Il26ll when studying the Glauber 
dynamics for the Ising model on regular trees. With this graph setting and in pres- 
ence of (+) -boundary condition, they proved in fact that the mixing time remains 
of order logn also at low temperatures (phase coexistence region), in contrast to 
the free boundary case where it grows polynomially in n IfTTl [4D . 

In this paper we extend the above result to a class of nonamenable graphs which 
includes trees, but also hyperbolic graphs with sufficiently high degree, and some 
suitable constructed expanders. Specifically, we consider the dynamics on an n- 
vertex ball of the graph with (+) -boundary condition, and prove that the spectral 
gap is Q(l) (i.e. bounded away from zero uniformly in ri) for all low enough tem- 
peratures and zero external field. This implies, by classical argument (see, e.g., 
H28ll ), an upper bound of order n on the mixing time. Notice that this result is in 
contrast with the behavior of the free boundary dynamics on hyperbolic graphs, for 
which the spectral gap is decreasing in n for all low temperatures, and bounded 
below by n^ a ^\ with exponent a((3) arbitrarily increasing with (3 IfTTl I4H. More- 
over, we give an example of an expander, in the above class of graphs, for which we 




(1.2) 
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prove that the mixing time of the free boundary dynamics is at least exponentially 
large in n. This provides a first rigorous example of graph where the mixing time 
shrinks from exponential to linear in n while passing from free to (+) -boundary 
condition. 

We remark that what we believe to be determinant for the result obtained in 
lf26ll for the dynamics on trees, is in fact the nonamenability of the graph. On 
the other hand, the possible presence of cycles, which are absent on trees, makes 
the structure of some nonamenable graphs more similar to classical lattices. Our 
results show that cycles are not an obstacle for proving the influence of boundary 
conditions on the mixing time. 

The work is organized as follows. In section 2 we give some basic definitions 
and state the main results. In section 3 we analyze the system at equilibrium and 
prove a mixing property of the plus phase. Then, in section 4, we deduce from this 
property a lower bound for the spectral gap of the dynamics and conclude the proof 
of our main result. Finally, in section 5, we give an example of a graph satisfying 
the hypothesis of the main theorem, and prove for it an exponential lower bound 
on the spectral gap for the free boundary dynamics. 

2. The model: definitions and main result 

2.1. Graph setting. Before describing the class of graphs in which we are inter- 
ested, let us fix some notation and recall a few definitions concerning the graph 
structure. 

Let G = (V, E) be a general (finite or infinite) graph, where V denotes the vertex 
set and E the edge set. We will always implicitly assume that G is connected. The 
graph distance between two vertices x, y G V is defined as the length of the shortest 
path from x to y and it is denoted by d (x, y). If x and y are at distance one, i.e. if 
they are neighbors, we write x ~ y. The set of neighbors of x is denoted by N x , 
and \N X \ is called the degree of x. 

For a given subset S C V, let E(S) be the set of all edges in E which have both 
their end vertices in S and define the induced subgraph on S by G(S) := (S, E(S)). 
When it creates no confusion, we will identify G(S) with its vertex set S. 
For S C V let us introduce the vertex boundary of S 

8 V S = {x eV\S : By G S s.t. x ~ y} 

and the edge boundary of S 

d E s = {e = (x, y)eE s.t. x G S ,y e V \ S} . 

If G = (V, E) is an infinite, locally finite graph, we can define the edge isoperimetric 
constant of G (also called Cheeger constant) by 

z e (G) := inf ; S C V finite j . (2.1) 

Definition 2.1. An infinite graph G = (V,E) is amenable if its edge isoperimetric 
constant is zero, i.e. if for every e > there is a finite set of vertices S such that 
\8eS\ < e\S\. Otherwise G is nonamenable. 

Roughly speaking, a nonamenable graph is such that the boundary of every sub- 
graph is of comparable size to its volume. A typical example of amenable graph is 
the lattice Z d , while one can easily show that regular trees, with branching number 
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bigger than two, are nonamenable. We emphasize that nonamenability seems to 
be strongly related to the qualitative behavior of models in statistical mechanics. 
See, e.g., lfT5l [T9l l20l [29ll for results concerning the Ising and the Potts models, 
and flU [71 [HI O for percolation and random cluster models. 

In this paper we focus on a class of nonamenable graphs, that we call growing 
graphs, defined as follows. Given an infinite graph G = (V, E) and a vertex o eV, 
let B r (o) denote the ball centered in o and with radius r G N with respect to the 
graph distance, namely the finite subgraph induced on {x G V : d (o, x) < r}, and 

let L r (o) := {x G V : d(x, o) = r} = dyB^o). 

Definition 2.2. An infinite graph G = (V, E) is growing with parameter g > and 
root o EV, if 

min {\N X n L r+1 (o)\ - \N X n B r {o)\} = g . (2.2) 

We call G a (g, o)-growing graph. 

It is easy to prove that a growing graph in the sense of Definition 12.21 is also 
nonamenable. The simplest example of growing graph with parameter g, is an 
infinite tree with minimal vertex degree equal to g + 2, where the growing property 
is satisfied for every choice of the root on the vertex set. On the other hand, there 
are many examples of growing graphs which are not cycle-free. Between them 
we mention hyperbolic graphs, that we will prove to be growing provided that the 
vertex-degree is sufficiently high. 

Hyperbolic graphs are a family of infinite planar graphs characterized by a cycle 
periodic structure. They can be briefly described as follows (for their detailed 
construction see, e.g., 112211 . or Section 2 of ref. H2"7lP . Consider a planar graph 
in which each vertex has the same degree, denoted by v, and each face (or tile) is 
equilateral with constant number of sides denoted by s. If the parameters v and 
s satisfy the relation (v — 2)(s — 2) > 4, then the graph can be embedded in the 
hyperbolic plane H 2 and it is called hyperbolic graph with parameters v and s. It 
will be denoted by M(v, s). The typical representation of hyperbolic graphs make 
use of the Poincare disc that is in bi-univocal correspondence with H 2 (see Fig. 

ED. 

Hyperbolic graphs are nonamenable, with edge isoperimetric constant explicitly 
computed in [fTTT] as a function of v and s. Moreover, the following holds: 

Lemma 2.3. For all couples (v, s) such that s > 4 and v > 5, or s = 3 and v > 9, 
M(v, s) is a (g, o)-growing graph for every vertex o G V and with parameter g = 

g(v,s). 

The proof of this Lemma is postponed to Section 5, where we will also construct 
a growing graph that will serve us as further example of influence of boundary 
conditions on the mixing time. 

Let us stress, that due to the possible presence of cycles in a growing graph, a 
careful analysis of the correlations between spins will be required. This is actually 
the main distinction between our proof and the similar work on trees Il26l] . 

2.2. Ising model on nonamenable graphs. The Ising model on nonamenable 
graphs has been investigated in many papers (see, e.g, If20ll for a survey). A general 
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FIGURE 2.1. The hyperbolic graph H(4, 5) in the Poincare disc representation. 

result, concerning the uniqueness/non-uniqueness phase transition of the model, 
is the following ITT5H : 

Theorem 2.4 (Jonasson and Steif). If G is a connected nonamenable graph with 
bounded degree, then there exists an inverse temperature j3 > 0, depending on the 
graph, such that for all (3 > f3 there exists an interval of h where G exhibits a phase 
transition. 

Thus, contrary to what happens on the Euclidian lattice, the Ising model on 
nonamenable graphs undergoes a phase transition also at non zero value of the 
external field. 

Though some properties of the Ising model are common to all nonamenable 
graphs, the particular behavior of the system may differ from one family to another 
one, also depending on other geometric parameters. Since we will be especially 
interested in hyperbolic graphs, we recall briefly the main results concerning the 
Ising model on these graphs, and stress which are the main differences from the 
model on classical lattices. 

It has been proved (see fl30], [33], [34]]) that the Ising model on H(t>, s) exhibits 
two different phase transitions appearing at inverse temperatures f3 c < f3' c . The 
first one, f3 c , corresponds to the occurrence of a uniqueness/non-uniqueness phase 
transition, while the second critical temperature refers to a change in the properties 
of the free boundary condition measure \J . Specifically, it is defined as 

p% := inf{/3 > p e : // = + /0/2} , (2.3) 

where u. + and u~ denote the extremal measures obtained by imposing, respectively, 
(+)- and (— )-boundary condition. As is well explained in H34l (see also fl8]] for 
more details), using the Fortuin-Kasteleyn representation it is possible to show 
that (3' c < oo for all hyperbolic graphs, in contrast to the behavior of the model on 
regular trees where pf ^ (p + + p~)/2 f° r an finite (3 > j3 c . From Definition 12.31 
it turns out that for (3 C < (3 < (3' c , when this interval is not empty (see U3410 , the 
measure pf is not a convex combination of and p~. This implies the existence 
of a translation invariant Gibbs state different from and p ~ , in contrast to what 
happens on Z d ||2fl. 
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Another interesting result concerning the Ising model on hyperbolic graphs, is 
due to Sinai and Series If30ll . For low enough temperatures and h = 0, they proved 
the existence of uncountably many mutually singular Gibbs states which they con- 
jectured to be extremal. This points out, once more, the difference between the 
system on hyperbolic graphs and on classical lattices, where it is known that the 
extremal measures are at most a countable number. 

In this paper we are interested in the region of the phase diagram where the 
dynamics is highly sensitive to the boundary condition, namely when the tempera- 
ture is low and the magnetic field is zero (phase coexistence region) . Let us explain 
the model in detail and give the necessary definitions and notation. 

Let G = (V, E) be an infinite (g, o)-growing graph with maximal degree A. For 
any r G N, we denote by B r = (V r , E r ) c G the ball with radius r centered in o. 
When it does not create confusion, we identify the subgraphs of G with their vertex 
sets. Given a finite ball B = B m and an Ising spin configuration r G Vl G , let Q T B c 
{±l} Bu5vB be the set of configurations that agree with r on d v B. Analogously, for 
any subset A C V m and any rj g Vt T B , we denote by Vt\ c {±i} Au9 v A the set of 
configurations that agree with 77 on dyA. The Ising model on A with 77-boundary 
condition (b.c.) and zero external field is thus specified by the Gibbs probability 
measure n A , with support on Q A , defined as 

AS?) = ^7m ex P(^ Yl a ^y)^ C2.4) 

(x,y)£E(A) 

where Z((3) is a normalizing constant and the sum runs over all pairs of nearest 
neighbors in the induced subgraph on A = A U dyA. 

Similarly, the Ising model on A with free boundary condition is specified by the 
Gibbs measure supported on the set of configurations VLa ■= {±1} A . This is 
defined as in A2.4D by replacing the sum over E(A) in a sum over E(A), namely 
cutting away the influence of the boundary dyA. Notice that when A = V m , \iy m 
is simply the Gibbs measure on B with boundary condition r (77 agrees with r on 
d v V m = d v B) and /j lVm is the Gibbs measure on B with free boundary condition. 

We denote by Ta the a-algebra generated by the set of projections {t[ x } x <^a from 
{±1} A to {±1}, where n x : a 1— ► a x , and write / e T A to indicate that / is JF A - 
measurable. Finally, we recall that if / : VT B — > R is a measurable function, the 
expectation of f w.r.t. fi v A is given by fi v A (f) = J^aen^Ai^fi 17 ) an( ^ tne variance 
of f w.r.t. fjP A is given by Var A = ^\(f 2 ) ~ ^\{f) 2 - ^ e usually think of them as 
functions of 77, that is fiA(f)(v) = ^\{f) an ^ ^^A{f){v) = Var^(/). In particular 
HaU) , Var j4 (/) G T A -- 

In the following discussion we will be concerned with the Ising model on B with 
(+)-b.c. and we will use the abbreviations T and fi instead of T-g and /j, b , 
and thus fi(f) and Var(/) instead of fj, B {f) and Var^(/). 

2.3. Glauber dynamics and mixing time. The Glauber dynamics on B with (+)- 
boundary condition is a continuous time Markov chain (cr(£)) t > on fi + with Markov 
generator C given by 



(£/)(*) = E^'M**) - > 

xeB 



(2.5) 
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where o x denotes the configuration obtained from a by flipping the spin at the 
site x and c x (a) is the jump rate from a to o x . We sometimes prefer the short 
notation V x f(cr) = [f(cr x ) — f(a)]. The jump rates are required to be of finite- 
range, uniformly positive, bounded, and they should satisfy the detailed balance 
condition w.r.t. the Gibbs measure /i. Although all our results apply to any choice 
of jump rates satisfying these hypothesis, for simplicity we will work with a specific 
choice called heat-bath dynamics: 

c x (a) := v°(a x ) = - 1 - where w x {&) := exp(2/3<r x Y]<r v ). (2.6) 

It is easy to check that the Glauber dynamics is ergodic and reversible w.r.t. the 
Gibbs measure \i, and so converges to p by the Perron-Frobenius Theorem. The 
key point is now to determine the rate of convergence of the dynamics. 

A useful tool to approach this problem is the spectral gap of the generator C, that 
can be defined as the inverse of the first nonzero eigenvalue of C. 

Remark 2.5. Notice that the generator C is a non-positive self-adjoint operator on 
£ 2 (f2 + ,/i). Its spectrum thus consists of discrete eigenvalues of finite multiplicity that 
can be arranged as = A > — Ai > — A 2 >...,> — Xn-i, if |^ + | = N, with Aj > 0. 

An equivalent definition of spectral gap is given through the so called Poincare 
inequality for the measure \i. For a function / : fi + i— > R, define the Dirichlet form 
of / associated to C by 

V U) ■■= ^E^( c -[ v ^] 2 ) = E^ Vax -(/))> (2 - 7) 

xGB xeB 

where the second equality holds under our specific choice of jump rates. The 
spectral gap of the generator, c gap (n), is then defined as the inverse of the best 
constant c in the Poincare inequality 

Vax(/) < cV(f) , Wf e f 2 (Q + , fi), (2.8) 

or equivalently 

c gap (fx) := inf {^^y; Var(/) ± o} . (2.9) 

Denoting by P t the Markov semigroup associated to C, with transition kernel 

Pt(cr, rj) = e tc (a, rj), it easy to show that 

Var(P t /) < e~ 2c! > ap(fl)t Var(f). (2.10) 

The last inequality shows that the spectral gap gives a measure of the exponential 
decay of the variance, and justifies the name relaxation time for the inverse of the 
spectral gap. 

Moreover, let h a t denote the density of the distribution at time t of the process 
starting at a w.r.t. p,, i.e. h? t (rj) = p ^vJ" . For 1 < p < oo and a function / e 
p,), let \\f\\p denote the £ p norm of / and define the time of convergence 

t p = min ji > : sup \\h° - l|| p < e _1 | , (2.11) 

that for p = 1 is called mixing time. A well known and useful result relating t p to 
the spectral gap (see, e.g., ||28l0 , when specializing to the Glauber dynamics yields 
the following: 
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Theorem 2.6. On an n-vertex ball B c G with (r)-boundary condition, 

Cgap{^Y l <n< Cgap^fl)' 1 X CU , (2.12) 

where p, = fi T B and cis a positive constant independent of n. □ 

We stress that a different choice of jump rates (here we considered the heat-bath 
dynamics) only affects the spectral gap by at most a constant factor. The bound 
stated in Theorem 12.61 is thus equivalent, apart for a multiplicative constant, for 
any choice of the Glauber dynamics. 

Before presenting our main result, we recall that the Glauber dynamics for the 
Ising model on regular trees and hyperbolic graphs has been recently investigated 
by Peres et al. [fTTl 01 . In particular, they consider the free boundary dynamics on 
a finite ball B c G, G hyperbolic graph or regular tree, and prove that at all tem- 
peratures, the inverse spectral gap (relaxation time) scales at most polynomially in 
the size of B, with exponent a(f3) | °o as j3 — > oo. Let us stress again that under 
the same conditions, the dynamics on a cube of size n in the (/-dimensional lattice, 
relaxes in a time exponentially large in the surface area n^ d ~ x ^ d . 

2.4. Main results. We are finally in position to state our main results. 

Theorem 2.7. Let G be an infinite (o, g)-growing graph with maximal degree A. 
Then, for all (3 1, the Glauber dynamics on the n-vertex ball B with (+)-boundary 
condition and zero external field has spectral gap 

As a corollary we obtain that, under the same hypothesis of the theorem above, 
the mixing time of the dynamics is bounded linearly in n (see Theorem l2.6D . 

This result, applied to hyperbolic graphs with sufficiently high degree, provides 
a convincing example of the influence of the boundary condition on the mixing 
time. Indeed, for all (3 3> 1, due to the fact that the free boundary measure on 
M(v,s) is a convex combination of n + and p~ (see section [2721) . it is not hard to 
prove that the spectral gap for an n-vertex ball in the hyperbolic graph with free 
boundary condition, is decreasing with n. Most likely it will be of order n~ a ^\ with 
a{(5) t oo as j3 — > oo, as in the lower bound given in |fT7i l4"H. The presence of (+)- 
boundary condition thus gives rise to a jump of the spectral gap, and consequently 
it speeds up the dynamics. 
Remarks. 

(i) We recall that on 1/ not much is known about the mixing time when (3 > f3 c , 
h = and the boundary condition is (+), though it has been conjectured that 
the it should be polynomial in n (see [9] and H5H). 

(ii) A result similar to Theorem \2.7\ has been obtained for the spectral gap, and 
thus for the mixing time, of the dynamics on a regular b-ary tree (see Il26l ). 
In particular it has been proved that while under free-boundary condition the 
mixing time on a tree of size n jumps from \ogn to n e ^ when passing a 
certain critical temperature, it remains of order log n at all temperatures and 
at all values of the magnetic field under (+)-boundary condition. However we 
stress that while trees do not have any cycle, growing graphs, and in general 
nonamenable graphs, can have many cycles, as well as the Euclidean lattices. 
The theorem above can thus be looked upon as an extension of this result to a 
class of graphs which in some respects are similar to Euclidean lattices. 



GLAUBER DYNAMICS ON NONAMENABLE GRAPHS: BOUNDARY CONDITIONS AND MIXING TIME 9 

(iii) At high enough temperatures (one phase region) the spectral gap of the dy- 
namics on a ball Be G, where G is an infinite graph with hounded degree, 
is for all boundary conditions, as can be proved by path coupling tech- 
niques Il32l1 . This suggests that the result of Theorem \2.7\ should hold for all 
temperatures, as for the dynamics on regular trees, and not only for (3 3> 1. 
At the moment, what happens in the intermediate region of temperature, still 
remains an open question. 

The following result provides a further example of influence of boundary condi- 
tions on the dynamics. 

Theorem 2.8. For all finite j 6 N, there exists an infinite (g , o)- growing graph G 
with bounded degree, such that, for all j3 ^> 1, the Glauber dynamics on the n-vertex 
ball B with free boundary condition and zero external field has spectral gap O (e~ 6n ), 
with9 = 6(g,p)>0. 

Combining this with Theorem 12.61 and Theorem |2.7[ we get a first rigorous ex- 
ample where the mixing time jumps abruptly from exponential to linear in n while 
passing from one boundary condition to another. 

We now proceed to sketch briefly the ideas and techniques used along the paper. 
The proof of our main result, Theorem |2.7[ is based on the variational definition 
of the spectral gap and it is aimed to show that the Gibbs measure relative to the 
system satisfies a Poincare inequality with constant c independent of the size of 
B. We will first analyze the equilibrium properties of the system conditioned on 
having (+) -boundary, and under this condition we will deduce a special kind of 
correlation decay between spins. The proof of this spatial mixing property rests on 
a disagreement argument and on a Peierls type argument. 

The second main step to prove Theorem 12.71 is deriving a Poincare inequality 
for the Gibbs measure from the obtained notion of spatial mixing. This will be 
achieved by first deducing, via coupling techniques, a like-Poincare inequality for 
the marginal Gibbs measure with support on suitable subsets, and then iterating 
the argument to recover the required estimate on the variance. 

The proof of Theorem |2.8| is given by the explicit construction of a growing graph 
with the property of remaining "expander" even when the boundary of a finite ball 
is erased, as in the free measure. Using a suitable test function, we will prove the 
stated exponentially small upper bound on the spectral gap. 

3. Mixing properties of the plus phase 

In this section we analyze the effect of the (+) -boundary condition on the equi- 
librium properties of the system. In particular, we prove that the Gibbs measure 
/i = /i^ satisfies a kind of spatial mixing property, i.e. a form of weak dependence 
between spins placed at distant sites. 

Before presenting the main result of this section, we need some more notation 
and definitions. Recall that for every integer i, we denoted by B { = (V^ Ej) the ball 
of radius % centered in o, and by B = B rn the ball of radius m such that \V m \ = n. 
Let us define the following objects: 

(i) the i-th level L t = {x G V : d(x, o) = i} = dyB^i, 

(ii) the vertex-set F { C B given by F { := {x G Bf_ x C\ B} ; 

(iii) the cr-algebra T% generated by the functions ir x for x G F£ = B^i. 
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We will be mainly concerned with the Gibbs distribution on F, with boundary con- 
dition rj G which we will shortly denote by fj% = yiP F . = \r] G analo- 
gously we will denote by Vaif the variance w.r.t. 

Notice that {-Fj}™^ 1 is a decreasing sequence of subsets such that V m = F d 
F x d . . . d F m+1 = 0, and in particular = fi>i{f), for all finite i, and 

fi m+ i(f) = f. The set of variables {/J>i(f)}i>o is a Martingale with respect to the 
filtration {J^}i> . 

For a given % g {0, . . . , m} and a given subset S c L i} we set U = F i+l U S 
and consider the Gibbs measure conditioned on the configuration outside U being 
t G which as usually will be denoted by u-Jj. 

We are now able to state the following: 

Proposition 3.1. Let G a (g,o)-growing graph with maximal degree A. Then there 
exists a constant 5 = 5(A) > such that, for every (3 > j^, every r G and every 
pair of vertices x G S C Li and y G Li \ S, i G {0, . . . , m}, 

= +) - = +)l < ce"^^ , (3.1) 

with f3' := 2g[3 — 5 > and for some constant c > 0. 

Let us briefly justify the above result. Since the boundary of B is proportional 
to its volume, the (+)-b.c. on B is strong enough to influence spins at arbitrary 
distance. In particular, as we will prove, the effect of the (+) -boundary on a given 
spin a x , weakens the influence on a x coming from other spins (placed in vertices 
arbitrary near to x) and gives rise to the decay correlation stated in Proposition 
13.11 Notice that the correlation decay increases with j3. 

The proof of Proposition 13. II is divided in two parts. First, we define a suitable 
event and show that the correlation between two spins is controlled by the proba- 
bility of this event. Then, in the second part, we estimate this probability using a 
Peierls type argument. Throughout the discussion c will denote a constant which is 
independent of \B\ = n, but may depend on the parameters A and g of the graph, 
and on (3. The particular value of c may change from line to line as the discussion 
progresses. 

3.1. Proof of Proposition 13.11 Let us consider two vertices x G S C L, L and 
y G Li \ S, such that d(x,y) = £, and a configuration r G Let t v,+ be the 
configuration that agrees with r in all sites but y and has a (+)-spin on y; define 
analogously r y ~ and denote by and the measures conditioned on having 
respectively r y > + - and -b.c. With this notation and from the obvious fact that 
the event {a : o x = +} is increasing, we get that 

\t%(°* = +)- = +)l = »u + (°* = +)- = +) • (3-2) 

In the rest of the proof we will focus on the correlation in the r.h.s. of (13 .2D . 

In order to introduce and have a better understanding of the ideas and tech- 
niques that we will use along the proof, we first consider the case £ — 1, which is 
simpler but with a similar structure to the general case £ > 1. 

3.1.1. Correlation decay: the case £ = 1. Assume that £ = 1, namely that x and y 
are neighbors. Denoting by ^ the measure with (— )-b.c. on U c = B,i \ {S} and 
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(+)-b.c. on b\B, we get 

l4} + {(?x = +)- lAf^x = +) = Hu~{(T x = -)- Hu + {°x = -) 

< i4}~( a x = -) 

< /i^(a x = -), (3.3) 

where the last inequality follows by monotonicity. The problem is thus reduced to 
estimate the probability of the event {a : o x = — } w.r.t. /x^. 

Let fC be the set of connected subsets of U containing x and write 

K = \_\ Kp with }C P = {C g fC s.t. \C\ = p} . 

p>i 

For any configuration a G we denote by the maximal negative component 
in fC admitted by a, i.e. 

With this notation the event {a : a x = — } can be expressed by means of disjoint 
events as 

{a: a x = -} = \_\ \_\ {a : = C} , (3.5) 

P >iceK p 

and then 

tot?* = -) = E E w = c ) • (3 - 6) 

Let us introduce the symbol o ~ C for a configuration cr such that a c = — and 
Cd v cr\u — +• The main step in the proof is to show the following claim: 

Claim 3.2. If G is a (g, o)-growing graph with maximal degree A, then, for any subset 

CcU, 

/V(<7 ~ C) < e" 29/3|c| . (3.7) 

The proof of Claim 13.21 is postponed to subsection 13.21 Let us assume for the 
moment its validity and complete the proof of the case £ = 1. By Claim [3721 and 
from the definition of K^ a \ we get 

H~(K^ =C)< e -^P\c\ . (3. 8) 

We now recall the following Lemma due to Kesten (see [TlolO . 

Lemma 3.3. Let G an infinite graph with maximum degree A and let C p be the set of 
connected sets with p vertices containing a fixed vertex v. Then \C P \ < (e(A + l)) p . 

Applying Lemma \33\ to the set fC p , we obtain the bound \fC p \ < e Sp , with 6 = 
1 + log(A + 1). Continuing from (13TB1) . we finally get that for all (3' = 2gf3 - 8 > 0, 
i.e. for all (3 > j-, 



= -) < E E 

p>l C&K P 



e -2gPp 



< Y^e~ 29(3p e 5p 
p>i 

< ce-P (3.9) 
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which concludes the proof of (13. ID in the case £ = 1 . 

Notice that the argument above only involves the spin at x, and thus applies for 
all pairs of x, y G Li, independently of their distance. Anyway when d(x, y) > 1 
this method does not provide the decay with the distance stated in Proposition 13. 11 
and a different approach is required. 

3.1.2. Correlation decay: the case £ > 1. Let us now consider two vertices x <E S C 
Li and y G Li \ S, such that d(x,y) = £ > 1. Before defining new objects, we 
want to clarify the main idea beyond the proof. Since the measure fixes the 
configuration on all sites in U c = Bi \ S, the vertex y can communicate with x only 
through paths going from x to y and crossing vertices in U. However, the effect of 
this communication can be very small compared to the information arriving to x 
from the (+)-boundary. In particular, if every path starting from y crosses a (+)- 
spin before arriving to x, then the communication between them is interrupted. 
Let us formalize this assertion. 

We denote by C the set of connected subsets C C U U {y} such that y G C, and 
call an element CeCa component of y. For every configuration a G we define 
as the maximal component of y which is negative on n U, i.e 

Observe that the spin on y is not fixed under the event {a : = C}. Finally, let 
C := {C g C s.t. x £ C} and define the event 

A := {a : C> } G C } = |J {a : C» = C} . (3.11) 

CeC» 

Then we have 

^-(a x = +\A) = J2h$-((t x = +,CW = C | A) 

E cgc ^r(^ = +,g (ff) = g) 
Ec ec ^r(c (<T) = ^) 

Ec e c^r(^ (<T) = c) 

> min = + I C (,t) = C) . (3.12) 

Notice that when the measure n^f is conditioned on the event {a : = C}, 
the spin configuration on dyC is completely determined by the boundary condi- 
tion: on dyC fi U it is given by all (+) -spins and on dyC n f/ c it corresponds to T y '~. 
Hence, spins on U \ (C U 9yC) become independent of spins on C, and we get 

$-(.\C^ = C) = fx y £(-\a z = +,zed v CnU) 

= ^ + (.\a z = +, ze (Cud v C)nU) 

> ^ + (-), (3.13) 

where the last inequality follows by stochastic domination. Being {a : a x — +} an 
increasing event, and from (I3.12D and (I3.13D . we get 

fi y u~(cr x = + \A)>ii y [ } + (a x = +), 
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which with the obvious fact that fif} (a x = +) > fif} (a x — + \ A) fxf} (A) , implies 

= +) - = +) < flr ( A °) ■ ( 3 - 14 ) 

By monotonicity and being A c a decreasing event, we get the inequality iif}~(A c ) < 
l^u(A c ) , where, we recall, [i^ denotes the measure on U conditioned on having all 
(— )-spins on U c . We now focus on /j^(A c ). 

Let C^ denote the set of components of y containing x, and for every pGN, let 
C p be the set of components in C^ with p vertices, i.e 

C p := {C E C+ % s.t. \C\ = p} := |J C p . 

p>0 

Notice that if C E C^ , then |C| > since d(x, y) = i. Thus, A c can be expressed 
by means of disjoint events as 

A ° = U U { a: c{a) = °}> (3 - 15) 

p>£+i ceCp 

and we get 

/%(^ c ) = E E ^ cia) = c ^ ■ (3 - 16) 

P >e+i cec p 

Since d v {C \ {y}) n U C <9yC n C/, we observe that the event {a : C (,j) = C} = {a : 

a c\{y} = - ' a d v cnu = +} is a subset of {a : crc\ fa} = - , cra v (C\{y})nu = +} = W '■ 
o ~ C \ {y}}. Applying the result stated in Claim 13721 to the set C \ {y}, we obtain 
the bound 

/i-(CW = CO < e-^^CIci- 1 ) , (3.17) 

which holds under the same hypothesis of the claim. Continuing from (I3.16D . we 
then have that for all j3' = 2g/3 - 5 > 0, i.e. for all (3 > ^, 

»uW < E E e " 9/3(p_1) 

p>e+i c&c p 
p>e 

< ce- 13 ' 1 , (3.18) 

where in the second line we used the bound \C P \ < e 5p due to Lemma [3T3l This 
concludes the proof of Proposition [370 In the next subsection we will go back and 
prove Claim [321 

3.2. Proof of Claim [3721 To estimate the probability /%(<r ~ C), we now appeal 
to a kind of Peierls argument that runs as follows (see also Ifl5l0 . Given a subset 
CCU, we consider the edge boundary OeC and define 

d+C := {e = (z,w) E d E C : z, w E U} 1Q1 
d-C := {e = (z,w) E d E C : z or w E U c } ' L J 

The meaning of this notation can be better understood if we consider a configura- 
tion a eVL^ such that = C (see (13.10D ). In this case a has (—) -spins on both 
the end-vertices of every edge in d-C and a (+)-spin in one end- vertex of every 
edge in d+C. Similarly if we consider a such that = C (see (I3.4D ). 
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For every a G fly such that cr ~ C, let a* G fi^ denote the configuration obtained 
by a global spin flip of a on the subset (7, and observe that the map a — > cr* is 
injective. This flipping changes the Hamiltonian contribute of the interactions just 
along the edges in d E C. In particular a* loses the positive contribute of the edges 
in d + C and gains the contribute of the edges in cLC, and then we get 

Hv(a*) = H v {a) - 2{\d + C\ - |cLC|) . (3.20) 

From this, we have 



fi v (a ~ C) = 



< 



_ e -2/3(|e+c|-|a_c|) ^ (3.21) 

where in the first inequality we reduced the partition function to a summation over 
{cr : cr ~ C} and then we applied (13.20D . 
The following Lemma concludes the proof of Claim I3T21 

Lemma 3.4. Let G a (g, o)-growing graph with maximal degree A. Then, for every 
subset CCU, 

\d+C\ - |<9_C| > g\C\ . (3.22) 

Proof. For a subset CCU, we define the downward boundary of C, c\C, as the 
edges of d E C such that the endpoint in C is in a higher level (strictly small index) 
than the endpoint not in C, i.e. 

dfi = {(u, v) e d E C : 3j s.t. u G Lj D C, v G . 

We then define the not-downward boundary of C, d?C, as the edges in d E C which 
are not-downward edges, i.e. d?C = d E C \ d^C. 

Notice that dyC C d+C, while d?C D d_C. In particular, inequality (13C221) follows 
from the bound 

|^t7| - |^t7| > ^|C| . (3.23) 
For all j > 0, define Cj = C n and notice that, by the growing property of G, 

I^C,-! - \d,C 5 \ > glCjl . (3.24) 
Moreover, one can easily realize that 

\ dl C\ - \dfC\ = \ 9 1 C j\ ~ l^il • ( 3 - 25 ) 
i>o 

In fact, all edges belonging to d E Cj, for some j, but not belonging to d E C, are 
summed once as downward edges and subtracted once as non-downward edges. 
In conclusion, the last two inequalities imply bound (I3.23D and conclude the proof 
of the lemma. □ 
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4. Fast mixing inside the plus phase 

In this section we will prove that the spectral gap of the Glauber dynamics, in 
the situation described by Theorem 12.71 is bounded from zero uniformly in the size 
of the system. From Definition 12.91 of spectral gap, this is equivalent to showing 
that for all inverse temperature j3 » 1, the Poincare inequality 

Var(/) < cV(f) , V/eL 2 (fl + ,7,,i) 

holds with constant c independent of the size of B. 

First, we give a brief sketch of the proof. The rest of the section is divided 
into two parts. In the first part, from the mixing property deduced in section 
[3] and by means of coupling techniques, we derive a Poincare inequality for some 
suitable marginal Gibbs measures. Then, in the second part, we will run a recursive 
argument that together with some estimates, also derived from Proposition 13.11 
will yield the Poincare inequality for the global Gibbs measure /j. 

4.1. Plan of the Proof. Let us first recall the following decomposition property of 
the variance which holds for all subsets D C C C B, 

Var^(J) = /4[Var D (/)] + Var^(/)] . (4.1) 

Applying recursively Q4.1D to subsets B = F D F 1 z> . . . D F m+ i = and recalling 
the relations ^(^ +1 (/)) = ^(/) and fi m +i(f) = f, we obtain 

Var(J) = /x[Var m (/)] + Var[/x m (/)] 

= //[Var m (/i m+ i(/))] + /i[Var m _i(/! m (/))] + Var[// m _i(// m (/))] 



= ^//[Vari^ (/))]. (4.2) 

Notice that A4.2D can also be seen as a decomposition of the Martingale given by 
the set of variables {yUj(/)}j>o respect to the filtration {^}j>o- 

To simplify the notation we define gi := j^i(f) for all z = 0, . . . , m + 1. Notice 
that gi G T{. Inserting in (I4.2D , we then have that 

m 

Var(/) = ^/x[Va ri (^ +1 )]. (4.3) 

i=0 

The proof of the Poincare inequality for fi, with constant independent of the size 
of the system, is given in the following two steps: 

(1) Proving that Vr 6 Q + and i e {0, . . . , m}, there exist suitable vertex-subsets 
{K x } x( z Li , K x 3 x, such that the like-Poincare inequality 

Var[(ft + i) < c ^(Var^(^ +1 )) (4.4) 

holds with constant c uniformly bounded in the size of Li, 

(2) Relating the variance of g^ = ^(f) to the variance of / in order to get an 
inequality of the kind 

m m 

J2 X) ^( Var ^(^+0) < cP(/) +eJ2Yl ^(Varj^foi+i)) (4.5) 

2=0 a;GLi i=0 a;GLi 
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with e a small quantity for f3 3> 1 . 
Notice that from (I4.5D the inequality 

in 

E M(Var^(^ +1 )) < c(l - ^(J) 

follows with c(l - e)- 1 = fi(l) for all /3 > 1. Together with Eqs. (1431) and (l4~4lh 
this will establish the required Poincare inequality for /i and therefore will conclude 
the proof of Theorem 12.71 

4.2. Step 1: From correlation decay to Poincare inequality. In this section we 
prove that under the same hypothesis of Proposition 13. 11 the marginal of the con- 
ditioned Gibbs measure on some suitable subsets, satisfies a Poincare inequality 
with constant independent of the size of these subsets. 

To state the result, let us fix a subset S C Li and a configuration r G f2 + . We 
then define the measure 

v T s(p) ■= E Kv\TeF Bi \s), (4.6) 

which is the marginal of the Gibbs measure Hp uS on S, and denote by Var^r the 
variance w.r.t. v T s . We state the following: 

Theorem 4.1. For all (3 » 1 and for every subset S C Lj, r G f2 + and / G 

L 2 {Vl, Fs, Vs), the measure v T s satisfies the Poincare inequality 

Var u r(f) <c J2 Var^/)) . (4.7) 
with c = co(P,A,g) = l + 0(e-«*»). 

Remark 4.2. Before proceeding with the proo f o f Theorem |4T1 we point out that this 
resuit includes, as a particular case, inequality (\4.4\) for subsets K x = F i+1 U {x}. 
To see that, choose S = L t so that //(• {Fb^s) = A 4 *- -An easy computation shows 
that, for every function f G J^+i, !/£.(/) = //[(/) and ^(Var^/)) = //[( Far^(/)). 
Inequality (\4.7& then corresponds to the like-Poincare inequality 

VarJ(g i+1 ) < c ^ /4( ^ a? x(^+i)) , 
which concludes the first step of the proof of Theorem \2.7\ 

4.2.1. Proof of Theorem \4. 1 1 The proof of Theorem [4J] rests on the so called cou- 
pling technique. This is a useful method to bound from above the mixing time of 
Markov processes, introduced for the first time in this setting by Aldous [1] and 
subsequently refined to the path coupling ||3ll2Tll . See also [18] for a wider discus- 
sion on the coupling method. 

A coupling of two measure p,\ and /j, 2 on is any joint distribution p on Vl x Q 
whose marginal are \i\ and \x 2 respectively. Here, we want to construct a coupling 
of two Glauber dynamics on Vt s with same reversible measure v T s but different 
initial configurations. We denote by Cs the generator of this dynamics. We also 
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recall that for all 77 G Q T S , x G S and a G {±1}, the jump rates of the heat bath 
version of the dynamics (see Def. 12.61 ) are given by 

c x (r], a) = v T s {a x = a\rj G Fs\x) 

= n{a x = a 1 77 G .Fs^, r G Fb^s) 

where in the second line we applied the definition of v T s and used that {a x = a} G 
Fs, and in the last line we adopt the notation K x := F i+X U {x}. 

We now consider the coupled process (r/(t), £(t))t> on f2 5 xf2s defined as follows. 
Given the initial configurations (77, £), we let the two dynamics evolve at the same 
time and update the configurations at the same vertex. We then chose the coupling 
jump rates c x ((r],a), (£,&)) to go from (77, £) to {rf ,a ,£, x,h ), with a,b G {±1}, as the 
optimal coupling (see HT8TI ) between the jump rates fi v Kx {<J x = a) and /i^ K {a x = b). 
More explicitly, for a G {±1}, they are given by 

f c x {{r], a), (f , a)) = min{//^ (a x = a) ; = a)} ( 4 Q ) 

I c x ((r/, a), (f , -a)) = max{0 ; = a) - ^ Kx (a x = a)} 

We denote by C the generator of the coupled process, and by P t the correspondent 
Markov semigroup. Notice that from our choice of coupling jump rates, we get that 
the probability of disagreement in x, after one update in x of (77, £), is given by 

pLM ■■= = +)- /4>* = +)l • (4.io) 

Let us now consider the subset F c fl s x O s given by all couples of configura- 
tions which differ by a single spin flip in some vertex of S. One can easily verify that 
the graph (fi 5 , H) is connected and that the induced graph distance between con- 
figurations (77, £) gO s x Vl s , D(r],^), just corresponds to their Hamming distance. 
Let us also denote by ~E v ^[D(rj(t),^(t))] the average distance at time t between two 
coupled configurations of the process starting at (77, £). We claim the following: 

Claim 4.3. For all (3 3> I, there exists a positive constant a = a(/3,g, A) such that, 
for every initial configurations (77, £) G H, the process (77(7;), £(;£)) t > satisfies the in- 
equality 

j t E^[D( V (t),m)}\t=o <-«. (4.11) 

Proo f o f Claim 14.31 The derivative in t of the average distance, computed for t = 0, 
can be written as 

j t ^[D( v (t),m]\t=o = j t (PtD)(rj,0\t=o = (£D)(r),0 

= E E c x ((v,a)^,b))[D(v x ' a ,e' b )-D( V ,0\- (4-12) 

x£S a,fee{±l} 

Since (77, £) G -ff, there exists a vertex y & S such that i = rf. If x = 7/, then 
Pdi S (Vi V x ) — an d the distance between the updated configurations decreases of 
one. While if x ^ y, with probability Pd is (r],ri y ) the updated configurations have 
different spin at x and their distance increases by one. Continuing from (14.12D we 
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get, for all /3 > 1, 

~dt 



16S 



e>i 

< -(l-ce-^), (4.13) 
where in the second line we used the bound 

Pl(v,V y ) = l/4>* = +) - A#>» = +)l < ce-^) , (4.14) 

which holds for all /3' = 5(3 — 2g > as stated in Proposition [37TJ Claim [431 follows 
taking a = (1 — ce _/3 ') and /3 sufficiently large. □ 

Using the path coupling technique (see [3]) we can extend the result of Claim 
|4.3| to arbitrary initial configurations (rj, £) G fig x Q s , and obtain 

j t E^[D( V (t)^(t))}\ t=0 <-aD( V ,Z). (4.15) 

From (147151) it now follows straightforwardly that E^[D (77(f) , f (i))] < e- at D(r], £), 
and then we get 

P(77(t) 7^ £(*)) < E^(£>(r7(*),e(t))) < e- at D{ v ,0 . (4.16) 

To bound the spectral gap c ga p{v T s ) of the dynamics on S, we consider an eigen- 
function / of £ s with eigenvalue —Cg ap {y T s ), so that 

Kf(v(t)) = e t£s f(rj) = e- c ^s)tf( v ). 

Since the identity function has eigenvalue zero, and therefore it is orthogonal to /, 
it follows that v T s (f) = and Vg(E^f(^(t))) = 0, where v T s is the invariant measure 
for £ 5 . From these considerations and inequality (I4.16D , we have 

< 2||/i| 0O su P p(^)^ew) 

< 2||/|U^| e - a *. (4.17) 

From the last computation, which holds for all 77 G Q T S and for all t, we finally 
obtain that c gap {vg) > a independently of the size of S, which implies the Poincare 
inequality (14. 7D with constant c = a' 1 = 1 + 0(e -c/3 ). This concludes the proof of 
Theorem 073 □ 

4.3. Step 2: Poincare inequality for the global Gibbs measure. With the previ- 
ous analysis we obtained a like-Poincare inequality for the marginal of the measure 
m on the level L { (see Remark [472]) , which inserted in formula (14.31) provides the 
bound 

m 

Var(/) < c E ^ ^(Var^( ft+1 ))l • (4-18) 

i=0 a;GLi 
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Using the same notation as in Il26)1 . let us denote the sum in the r.h.s. of (14.181) 
by Pvar(/)- The aim of the following analysis is to study Pvar(/) i n order to find 
an inequality of the kind Pyar(/) < cD(f) + ePvar(/), with e = e((3,g,A) < 1 
independent of the size of the system. This would imply that 

Var(/) < c • P V ar(/) < <*Yrf>(f) - 

and then would conclude the proof of Theorem 12.71 

In this last part of the section, we will first relate the local variance of $ = 
with the local variance of /. This will produce a covariance term that will be 
analyzed using a recursive argument. 

4.3.1. Reduction to covariance. In order to reconstruct the Dirichlet form of / from 
Q4.18D , we want to extract the local variance of / from the local variance of g i+1 . 
Notice that w.r.t. the measure jj, Kx , the function g i+ i just depends on x. Fixing 

x e Li and r e Q + , and defining p{r) : = fi T Kx ((T x = +) and q{r) := H T Kx {cr x = -), 
we can write 



to {Vai Kx {g i+ i)) = ^2^i{r)p(T)q{T){V x g i+1 (T)) 2 . 



(4.19) 



Using the martingale property g i+ \ = ii i+ x{g i+2 ), the variance Var^ (ft+i) can be 
split in two terms, stressing the dependence on x of g i+2 and of the conditioned 
measure Let us formalize this idea. 

For a given configuration refl + we introduce the symbols 

if // / x 



T 



r, 



if y 



x 



T 



ify^x 
if y = x 



and define the density 

hJa) 



t+i 



cr 



with t4 +1 (h x ) = 1 



Whit this notation and continuing from (14.19D , we get 

Hi(V&r Kx (g i+1 )) = ^^i(T)p(T)q(T)[V x ^i + i(gi +2 )(T)} 2 

T 
T 

= ^2vi(r)p(T)q(T) //]"+! (V x pi+ 2 ) - vl+xihx, g i+2 ) 

T 

< 2j2Mp(r)q(r) 



/^[+i(Vxft+2) + [v T i+i(h x ,g i+ 2 



(4.20) 



(4.21) 



Consider the first term of A4.21D and notice that x gi+ 2 ) — ttf+ifflxf)- To 

understand this fact, it is enough to observe that the dependence on x of g i+2 = 
Hi +2 (f) comes only from /, since the b.c. on B i+ i is fixed equal to r + . Replacing 
^7 x gi+2 by V x f and applying the Jensen inequality, we get 
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Now, notice that p(r)^ 1 (tr + ) = = n%(a x = +) (n r K Ja + ) + ^(a )). 

This, together with the fact that V x / does not depend on x, means that the expres- 
sion on the r.h.s. of (14.221) equals to 

W (r)g(T) S/iL(ffK(^ = +) (VJWf < X) WW mf {—'' 



Since a agrees with r on K x , then 

q(t) VkSp* = -) VkM) 



K(°x = -) V t kS? x = -|tr) /^H "* 

g(r) 

;<t 3: = 

To bound l|/iJ|oc» we first write 



Ox = - 

(4.23) 



< h x (a) , (4.24) 



and hence, inf g <! } < HMoo 



Taking the supremum over a oih x , we then have 

IIMoo < 7< ^ 1 _ . r, (4.26) 

= +; e N x n L i+1 ) l - L^nz^ = -) 

where we denoted by the measure conditioned on having all minus spins in 
Bi and plus spins in d v B. Inequality (I3.9D , applied to U = F i+ i, implies that 

H*i(Pv = -)<<*-* , (4.27) 
with (3' = 2g/3 - 5 as in Proposition [370 Combining (14.261) and (l4~27l) , we get that 
for all /?>^, 

||M~ < 1 + ce"^'. (4.28) 
Altogether, inequalities fl4.22D - fl4.28D imply that 

J2^(r)p(r)q(r) (W+i(V^ +2 )) 2 < c lfH (Var^/)) , (4.29) 

T 

with ci = ci(/3, A, g) = l + 0(e _c/3 ). Thus, summing both sides of (I4.21D over x e L { 
and i e {0, . . . , m}, and applying inequality Q4.29D , we obtain 



Pvar(/) < 2 Cl £>(/) + 2 



i=0 xgLi 



^^(r)p(r)g(r) U T i+1 (h x , g i+2 



(4.30) 



Notice that since g m +2 = / is constant w.r.t. /i m+ i, then ^ +1 (h x , g m+2 ) = and 
the value m can be removed from the summation over i in the r.h.s. of (I4.30D . 
It now remains to analyze the covariance /J,l+x(h X} g^)- 

4.3.2. Recursive argument. Before going on with the proof, we need some more 
definitions and notation. For every x e L i} let D x denote the set of nearest neigh- 
bors of x in the level L i+1 (descendants of x). Given x e L { and i e N, let us define 
the following objects: 

(i) D x i :={?/€ L i+ i : d(y, D x ) < £} is the ^-neighborhood of D x in 

(ii) JF^ := a [oy : ye B i+i \ D x j) is the cr-algebra generated by the spins on 

Bi+i \ D x f, 

(iii) fi Xj t( ■ ) := jj (■ | F x j) is the Gibbs measure conditioned on the a-algebra T x $. 
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We remark that D xfi = D x , and that there exists some £ < 2(i + 1) such that, for 
all integers i > £ , D x/ = L i+1 and /i x/ = fi i+1 . 

We also remark that for any function / e L 1 (fi, jF i+1 , /x), the set of variables 
{A i x,^(/)}^eN is a Martingale with respect to the filtration {Fx/}^,!,...,^- 



Let us now come back to our proof and recall the following property of the 
covariance. For all subsets D C C C B, 



(4.31) 



Since the support of strictly contains the support of yu x . , we can apply the 
property (I4.31D to the square covariance (n1 +1 (h x , g i+2 )) 2 appearing in (14.30D , in 
order to get 

(n T w (h x ,g i+2 )) 2 < 2(^ 1 (fi x ^h X} g i+2 ))) 2 + 2(^ 1 (^ x ^h x ),fi xfi (g i+2 ))) 2 . (4.32) 
The first term in the r.h.s. of (I4.32D can be bounded, by the Schwartz inequality, as 



(^+i(^x,o( h x, 9i+2))) 2 < nl +1 {Vax X;0 (h x )) ■ fi T i+1 {V&r xfi (g i+2 )) ■ 
The second term can be rearranged and bounded as follows: 



(4.33) 



lv> T i+ i(Vx,o(h x ), fx xfi (g i+2 ))Y 



i+l 



^xfl{h x ) — Hi+i(h x ), g i+2 ) 



ti+i (^^-i(^) ~ Hx,e(h x )) ,g i+2 



J=i 



- \L I i+i(l J 'x,e-i(h x ) ~ HxA h x)i9i+2) 



Uh-i {Px,i{Hx,i-\{h x ),gi+ 2 )) 



(4.34) 



where in the second line, due to the fact that \i x ^ = for some £ , we substi- 
tuted /J, x , {h x )-Hi + i(h x ) by the telescopic sum Eti(^,<-i(U:^,<(^)) • A PPly in g 
again the Cauchy-Schwartz inequality to the last term in (I4.34D , we get 



[mI+iC^o^x), Hx,o{gi+z) )] 2 < 



t=i 

To conclude the estimate on the covariance, it remains to analyze the three quan- 
tities appearing in (I4.33D and (14.35K : 

(i) //T'j (Wzi x/ (g i+2 )), for all £ = 0, 1, ... , £ ; 

(ii) ^[^(Var^o^a,)); 

(iii) K+i (y &I x,eX^x,i-i{h x ))), for all £ — 1, . . . , £ Q , 
We proceed estimating separately these three terms. 

First term: Poincare inequality for the marginal measure on D x ^. 

Let us consider the variance Va,r X: £(g i+2 ) appearing in (i). By definition, the function 



(4.35) 
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g i+2 depends on the spin configuration on B i+1 . Since the measure ^jp xl fixes the 
configuration on B i+ i \ D x ^, it follows that 

Thus, for every configuration rj G we can apply the Poincare inequality stated 
in Theorem (14. ID to Var^ e (g i+2 ), and obtain the inequality 

lA+i (Var^(# i+2 )) < c Mlfi( Var ^(^+2))) (4-36) 

with c = 1 + 0(e~ cf3 ) independent of the size of system. 
Second term: computation of the variance of h x . 

Notice that, from definition (14.20D , it turns out that h x only depends on the spin 
configuration on D x . In particular, for all rj which agrees with t~ on //j, (h x ) = 
and 

Vaxl (h x )<\\h x \\l-1. (4.37) 
Together with inequality G4.28D , this yields the bound: 

Vai%(h x ) < ce- 13 ' =: kp . (4.38) 
Third term: the variance of fi x> e-i(h x ). 

We now consider the variance Var^ e ([i x ,e-i(h x )), with rj G f2 + and £ > 1. Applying 
the result of Theorem 14.11 we obtain 

Vax x/ (n x> z-i(h x )) < c Y ^(Varx z (/i^-i(^))) 

= c Y f^li(y^K z {^-i(h x ))) , (4.39) 

2e-Di,A-Di,<-i 

where in the last line we used that the function /J, x ,e-i(h x ) does not depend on the 
spin configuration on Z^-i- 

Let z G D x ^ \ D X)i _i, and for any configuration ( G Vl r ' D , let us denote by ( + and 
(~ the configurations that agree with ( in all sites but z, and have respectively a 
(+)-spin and a (— )-spin on z. The summand in (14. 391) can be trivially bounded as 

^ x ,e(y^K z (^x,£-i(h x ))) < I sup (mC-iC 1 *) -^L-iO*)) 2 - (4-40) 

Moreover, by stochastic domination, the inequality ^ x ^i{h x ) > ^i^i(h x ) holds. 

Let z/(cr, a') denote a monotone coupling with marginal measures n^t-x an d /4 
We then have 

fj£e-i( h x) - /4,£-i W = Y v (°> a "> ~~ MO) 

< AH^Iloo max (v(a y = +) - v(o' y = +)) 

= Ay/ixlloo max^.iK = +) - = +))> (4.41) 
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where we used that the function h x only depends on the spins on D x . 

From Proposition [3TT1 with U = F i+2 U D x j, and since d(z, y) > d(z } D x ) > £, the 
probability of disagreement appearing in (14.411) can be bounded as 

tC-ifrv = +)- t£t-i(°v = +)< ce_/3 '' • H.42) 

Putting together formulas Q4.39D - Q4.42D , and applying inequality (14.281) . we ob- 
tain that for all rj e 

Varl^^iK)) < Kpe-W (4.43) 

with k' = c(l + 0(e" c/3 )) . 



Conclusion. 

Let us go back to inequalities Q4.33D and 04.351) ■ Applying bounds fl4.36DJ4.38D 
and (I4.43D , we get respectively 

yeD x 

• [^+i{^,o{h x ),^ (g i+2 ))} 2 < k'^i 2 e- 2p '" ^ /^(Var^Q^)), 

t=\ yeD x , e 

where we included in kp and k'^ all constants non depending on (3. 

For all f3 1, there exists a constant e = e(/3, g, A) = 0(e~ c/3 ) such that kp < e 
and k'^e' 13 ' 1 < k'^e' 13 ' < e. Substituting e in the inequalities above and summing 
the two terms as in (I4.32D , we obtain 

2 io 

jil +1 (h x , g i+2 ) J <£^e^' £ ^I+i( Var ^(^+2)) • 
£=0 jefl,,, 

Inserting this result in the second term of formula (14.30D and rearranging the sum- 
mation, we get 

m—1 

xi yi+2, 

i=0 xdLi 



m—1 



^ "J2J2J2J2 e^M^ Ky (g i+2 )) 

i=0 x£Li £=0 y£D x ^ 
m—1 £q 

< zJ2 E ^(Var^( ft+2 ))£e-^nW 

i=0 yGLi+i £=0 



(4.44) 



where in the last line we denoted by n(£) the factor which bounds the number of 
vertices x such that a fixed vertex y belongs to D Xj £. Since n(£) grows at most like 
A e , the product e~P' l n{l) decays exponentially with £ for all f3 ^> 1. Thus the sum 
over £ e {0, . . . , £ } can De bounded by a finite constant c which will be included in 
the factor e in front of the summations. Continuing from Q4.44D , we get 



m—1 



E5> 

i=0 x£Li 



xi g 



n+2, 



< ^E ^2^ Vai K y {gi+i)) 

i=l yeLi 

< ePvar(/). (4-45) 
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Inserting this result in H-30D and noticing that e = 0(e cf3 ) < 1 for /3 large enough, 
we obtain 

Pvar(/) < 2c{D{f) + eP V ar(/) =► Pvar(J) < ^ V(f) . 
Together with inequality (I4.18D , this implies that 

Var(/) < c P V ar(/) < c 2 V(f) , 

that is the desired Poincare inequality with c 2 = c 2 (/3, g, A) = 2(1 + 0(e~ c/3 )) inde- 
pendent of the size of the system . Notice also that the lower bound on the spectral 
gap, l/c 2 , increases with j3 and converges to 1/2 when (3 | oo. This concludes the 
proof of Theorem 12. 71 □ 

5. Influence of boundary conditions on the mixing time 

In this section we discuss two examples of influence of the boundary condition 
on the mixing time derived from Theorem |2.7i In particular, we first prove Lemma 
12.31 which implies the applicability of Theorem 12.71 to hyperbolic graphs with suf- 
ficiently high degree. Then we prove Theorem 12.81 providing an explicit example 
of growing graph which exhibits the behavior stated in the theorem. 

5.1. Hyperbolic graphs. In this subsection we want to prove Lemma [231 Before 
starting the proof, let us recall the following definition: 

Definition 5.1. The number of ends, 8{G), of a graph G = (V, E), is defined as 
£{G) := sup { number of infinite connected components of G\K} , 

KCV 

infinite 

where G\K denotes the graph obtained from G by removing the vertices which belong 
to K and the edges incident to these vertices. 

It is well known that hyperbolic graphs have one-end, as well as all the lattices 
Z d with d > 2. This property, together with the planarity and the cycle-periodic 
structure of the hyperbolic graphs, will be the main ingredient of the next proof. 

Proof of Lemma \2^3\ Consider a planar embedding of H(t>, s) and assume, without 
loss of generality, that the vertices of any ball B r around a fixed vertex o e V, are 
in the infinite face of the graph. 

With the same notation introduced in the previous sections, let us consider x G 
L r = d v B r _i. Let D x = N x n L r+1 be the set of "descendants of x", and define 

S x := N x n L r and P x := N x n L r _ x , 

so that N x = D x U S x U P x . Thus, when specialized to hyperbolic graphs M(v, s), 
the definition 12.21 of growing graph corresponds to the condition: 

min {D x -S x -P x }= min {v - 2(S X + P x )} = g . (5.1) 

xeL r ,reN xgL r ,reN 

Let us start verifying the following general properties: 

(i) Every vertex x G L r , has at least one neighbor in L r+1 , i.e. \D X \ > 1; 

(ii) Every vertex x G L r , has at most two neighbor in L r , i.e. \S X \ < 2; 

(iii) Every vertex x G L r , has at most two neighbors in L r _i, i.e. \P X \ < 2. 
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Property (i) is proved by contradiction. Assume the existence of a vertex x e L r 
such that D x — $. By the periodic structure of the graph it follows that there are at 
least v vertices in L r with no descendants, spaced out by vertices in L r connected 
to M(v, s) \ B r by at least one edge. Then, it exists a finite path, from x to another 
vertex in L r with no descendants, dividing two infinite components. This is in 
contradiction we the property of having one end, thus we conclude that \D X \ > 1. 

Property (ii) is also proved contradiction. Assume the existence of a vertex x e 
L r linked to three vertices in L r . Then there exist two faces, included in B r , both 
passing through x and one of its neighbors on L r , say y (see Fig. 15.21 left frame). 
Consequently, or D y ^ 0, contradicting the planarity of the graph or the property of 
having one-end, or D y = 0, contradicting property (i). We conclude that \S X \ < 2. 

Property (iii) is clearly satisfied if v — 3, as a consequence of property (i). Let 
v > 4 and proceed by induction. For all vertices x in the first level, we have 
obviously that \P X \ = 1. Assume that property (iii) holds for all vertices of L r _i, 
and consider a vertex x G L r . If we assume by contradiction that P x = 3, then there 
exist two faces, included in B r , both passing through x and one of its parents, say 
y (see Fig. 15.21 right frame). Consequently, the vertex y can not have neighbors 
in L r other than x, because this would contradict the planarity of the graph or 
the property of having one-end. On the other hand, if s > 3, y can not be linked 
to the other parents of x because this would create a triangular face. Globally, 
and by the inductive hypothesis that \P y \ < 2, it follows that \N y \ < 3, which is 
a contradiction. Then P x < 2, which completes the induction step for s > 3. If 
s = 3, y is connected with the two other parents of x to create two triangular 
faces. Thus, by the inductive hypothesis that \P y \ < 2, it follows that |iV y | = 5, 
which is a contradiction since v > 7 for all triangular tilings (recall the condition 
(v — 2)(s — 2) > 4). Then \P X \ < 2, and by induction this concludes the proof of 
property (iii) also for s = 3. 




FIGURE 5.2. The presence of more than two neighbors of x in the same 
level of x (left frame) or in the level above x (right frame) is not consistent 
with the structure of hyperbolic graphs. The red lines correspond to edges 
which are forbidden by the geometric properties of the graph. 



Properties (ii) and (iii) imply that S x + P x < 4. From (I5.1D . we then get that for 
all v > 9, H(t> , s) is a growing graph with parameter g > v — 8. Moreover, it can 
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be easily verified that if s = 3 (triangular tilings), then g — v — 8, since in this case 
S x + P x — 4 for some x e V . 

If s > 4, the above condition can be improved by proving the following proper- 
ties: 

(iv) If s — 2k, with k > 2, then for every x e L r , S x = 0; 

(v) If s = 2A; + 1, with k > 2, then for every x e L r , S x + P x < 2. 

Let us prove property (iv) by induction. For all vertices x in the first level, we 
have obviously that \S X \ =0, otherwise we would have a triangular face. Assume 
that property (iv) holds for all vertices of L r _i, and consider a vertex x 6 L r . If we 
assume that x has a neighbor in L r , say y, then there exists a face F, included in 
B r , passing through x, y, and the respective parents (see Fig. 15.31 left frame) . The 
remaining sides of F are all included in B r _i, and by the inductive hypothesis they 
can only connect vertices on subsequent levels. In particular, this face can only 
have an odd number of sides, which is in contradiction with the hypothesis that 
s = 2k. We conclude that S x = for all x e V. 

Property (v) is clearly satisfied if v = 3, as a consequence of property (i). Thus 
let v > 4 and assume by contradiction that + |Pa.| = 3, which means, coherently 
to properties (ii) and (iii), that S x = 2 and P x = 1 or viceversa. In both cases, there 
exist two faces, included in B r , passing through x and one of its parents, say y (see 
Fig. 15.31 central and right frames). Then the vertex y can not have neighbors in L r 
other than x, because this would create a triangular face, or would contradict the 
planarity of the graph or the property of having one-end. Analogously, y can not 
have neighbors in L r _ l5 because this would create a triangular or a square face, 
while s = 2k + 1 by hypothesis. Then we get that |A^| = 1 + \P y \ < 3, which is in 
contradiction with the hypothesis that v > 4. We conclude that \S X \ + \P X \ < 2. 




FIGURE 5.3. Forbidden patterns. In the left frame it is shown that S x / 
only if the number of side in each face is odd. The right and central frames 
show that the condition IS^I + \P X \ > 3 is possible only if s = 3. 

All together, properties (iii)-(v) imply that if s > 4, then S x + P x < 2. Thus, 
for all couples (v, s) such that s > 4 and v > 5, M(v, s) is a growing graph with 
parameter g < v — 4. It can be easily verified that if s > 4 then g — v — 4, since 
S x + P x = 2 for some x eV. □ 

Together with Theorem \2.7[ the above lemma implies that for all (v, s) such that 
M(v, s) is a growing graph and for all /3 > 1, the dynamics on an n-vertex ball of 
W(v, s) with (+)-boundary condition, has spectral gap uniformly positive in n and 
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mixing time at most linear in n. Comparing this result with the behavior of the 
dynamics with free boundary condition 0, we get a convincible example of the 
influence of boundary conditions on the mixing time. 

5.2. Expanders. In this subsection we prove Theorem l2.8l by providing an explicit 
example of growing graph which exhibits the behavior stated in the theorem. 

5.2.1. Construction. Let us start with some definitions. 

Let G = (V, E) be a finite graph. The edge isoperimetric constant of G, defined 
in (12. ID for infinite graphs, is given by 



For c > and k, n e N, we then have the following definition: 

Definition 5.2. A finite graph G = (V, E) is called an (n, k, c)-expander if it is regular 
with degree k, \V\ = n, and i e (G) = c. 

Notice that, since G is finite, i e (G) > if and only if G is connected. In particular, 
every connected fc-regular graph on n vertices is an expander for some c > 0. 
However, one is usually interested in a family of expanders, that is a sequence 
of (n, k, c)-expander graphs such that k is fixed, c > e > 0, and n — > oo. That 
c > e > ensures that these graphs are highly connected, or, in other terms, that 
the sequence of expanders converges to an infinite regular nonamenable graph. It 
is easy to see that this happens only if k > 3. On the other hand, for all k > 3, 
a sequence of (n, k, c) -expanders exists, as it has been proved by showing that k- 
regular random graphs on n vertices are expanders with probability tending to 1 as 
n — > oo (probabilistic method) . See lfT2ll for a nice survey on expanders and their 
applications. 

Let T A denote an infinite rooted tree with constant degree A and root o <eV, and 
assume that A > 6. Fix an integer d < A — 2, and connect the vertices on each level 
L r = {x G V : d{x, o) = r} of the tree in such a way that induced subgraph on L r 
is an (A(A — l) r_1 , k r , c r )-expander, with 3 < k r < d and isoperimetric constant c r 
uniformly positive in r . We denote the infinite graph obtained with this procedure 
byX A ' d (see Fig.[5T4l). 



FIGURE 5.4. Construction of X 8 ' 3 to first level. The red lines correspond 
to the edges of the expander in L\. 

Notice that X A d has degree D > A + d, and by construction it is (g, o)-growing 
graph with g > A — 2 — d. Moreover, one can easily realize that the isoperimetric 




(5.2) 
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constant i e (B r ) of the ball B r c X A ' d centered in o, is uniformly positive in r > 1, 
namely, there is a constant e = e(A, d) > such that i e (B r ) > e for all r e N. 



5.2.2. Free boundary dynamics finite balls of X A ' d . Let us consider the Glauber 
dynamics on the n-vertex ball B = B m of the graph X A ' d with free boundary 
condition. We use the same notation of section 12.31 but here ji denotes the Gibbs 
measure over Vl = Vl B = {±.1} B with free boundary condition, that is, for all cr e O, 

[ ■ ' (x,y)eE(B) 

With this notation, we can state the following result. 

Proposition 5.3. Let B the n-vertex ball of the graph X A ' d . Then, for all (3 » 1, 
the Glauber dynamics on B with free boundary condition and zero external field, has 
spectral gap 0(e~ 9n ), with 9 = 9(A,d,{3) > 0. 

Proof. For every o e Vt, let m B (a) := J2 x eB °x De tne magnetization of a configura- 
tion o, and define the following characteristic function: 

Wo}(<t) = | ifms((T) < • (5-4) 

Without loss of generality we can restrict the analysis to odd values of n. Since in 
this case m B (cr) ^ for all a, by symmetry it follows that Var(l{ ms>0 }) = 1/4 and 
then, from d2T9D , 

<W(/i) < 4:P (l{m s >0}) = 2 ^yU (Cj. [V x (l{ mfl>0 })] 2 ) . (5.5) 

Notice that Va;(l{ ms> o})(o") ^ only if a is such that \m B (<j)\ = 1, namely only 
if there are (n + 2)/2 spins in cr with the same value, and only if x is one of the 
(n + 2)/2 vertices with spin of the same sign of m B (a). Then we get 

V(t {mB>0} ) < ^ ni\m B \ = 1) = 2±2 /i(m fl = 1) . (5.6) 

For every configuration cr G f2, we define the subsets 

A + (a) := {x e B s.t. cr^ = +} and A~(a) = {x E B s.t. a x = -} . 

Notice that the condition m B (a) = 1 is satisfied if and only if | A + (cr) | = (n + l)/2 
and |A~(cr)| = (n — l)/2. Thus, any configuration a with magnetization equal to 
1, is univocally correspondent to a partition of the vertex set of B into two subsets 
S and T of size (n + l)/2 and (n — l)/2, respectively, such that = S and 

A _ (cr) = T. Let "P denote the set of these partitions, and define E(S, T) as the set 
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of edges between S and T. Then we have 

fi{m B = 1) = 



,y)£E(B) a x a y) 



< 



Z{(3) 

V eMP(\E(S)\ + \E(T)\-\E(S,T)\)) 
^ Z{0) 

v exp(P(\E(S)\ + \E(T)\-\E(S,T)\)) 
{s j^ p ex P ((3(\E(S)\ + \E(T)\ + \E(S,T)\)) 

= exp(-2p\E(S,T)\) . (5.7) 

{s,T)ev 

Notice that \E(S, T)| = \8eT\. In particular, since B has uniformly positive isoperi- 
metric constant, namely i e (B) > e > for all n G N, 

\E{S,T)\>e\T\ = *(n-l). 

Continuing from (15. 7D , we get 



n(m B = 1) < 2^ e 



0e(n-l) 
(S,T)EV 



n-1 J e 
2 / 

= ce -(^-l°g2)n ; (5. 8 ) 

where in the third line we approximated the binomial factor using the Stirling 
formula. Inserting this bound in (15. 6D we get that, for all (3 > there exists a 
positive constant 6 = 9 (A, d, (3) such that 

2>(l {ms >o } ) < ^ce~^- lo ^ n < e- 6n , (5.9) 

which implies, by (IB3D . that c flop (/i) = 0(e~ 9n ). □ 

Theorem 12.81 follows from the above proposition, since by construction the in- 
finite graph X A,d is growing with parameter g > A — d — 2. Putting together 
Theorems 12.71 and Proposition 15.31 we get the first rigorous example in which the 
mixing time is at most linear in n for the (+) -boundary condition, while at least 
exponential in n for the free boundary condition. 
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